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Hourglass-shaped magnetic excitation spectra have been detected in a variety of doped transition- 
metal oxides with stripe-like charge order. Compared to the predictions of spin- wave theory for 
perfect stripes, these spectra display a different intensity distribution and anomalous broadening. 
Here we show, based on a comprehensive modelling for La5/3Sri/3Co04, how quenched disorder in 
the charge sector causes frustration, and consequently cluster-glass behavior at low temperatures, 
in the spin sector. This spin-glass physics, which is insensitive to the detailed nature of the charge 
disorder, but sensitive to the relative strength of the magnetic inter-stripe coupling, ultimately 
determines the distribution of magnetic spectral weight: The excitation spectrum, calculated using 
spin waves in finite disordered systems, is found to match in detail the observed hourglass spectrum. 



Magnetic excitations, as measured by inelastic neutron 
scattering (INS), serve as a fingerprint of the underlying 
magnetic state. In layered transition-metal oxides, so- 
called "hourglass" spectra are frequently observed [1-10], 
consisting of intensity at incommensurate (IC) wavevec- 
tors at low energy which, with increasing energy, first 
disperses towards the antiferromagnetic (AF) wavevector 
(7r,7r), leading to a large-intensity "resonance", and then 
disperses outwards again. In the case of cuprate super- 
conductors, explanations both in terms of weak-coupling 
excitations of a Fermi liquid and of strong-coupling col- 
lective modes of an underlying stripe state have been 
proposed [11-15]. 

More recently, hourglass spectra have been detected in 
doped cobalt at es and manganites [9, 10]. Because these 
are insulators, the magnetic excitations can be conclu- 
sively linked to collective local-moment physics. More- 
over, charge order is known to occur, i.e., magnetic and 
non-magnetic ions form a pattern of unidirectional two- 
dimensional (2D) stripes [16]. Indeed, assuming mag- 
netic moments to be arranged in stripes and coupled 
by superexchange interactions allows one to qualitatively 
understand some gross features of the measured spectra. 

However, the experimental data display pronounced 
deviations from the spin- wave spectra of perfect stripes: 
(i) At low energies, outward dispersing excitation 
branches are essentially completely missing, (ii) The data 
typically show a large and anisotropic intrinsic broaden- 
ing, (iii) For a sizeable range of intermediate energies a 
"vertical dispersion" is observed near (7r,7r). (iv) Com- 
pared to perfect-stripe spin- wave results, spectral weight 
is redistributed from higher energies into the broad "res- 
onance" (or saddle-point) feature. It has been suggested 
that quenched disorder plays a role in explaining these 
features, but concrete calculations have been restricted 
to ad- hoc broadening of perfect-stripe results, and deeper 
insights are lacking [9, 10, 16]. 

In this paper we analyse spin- wave spectra of stripes in 
2D insulators, treating disorder exactly in large finite-size 
systems. Importantly, stripe disorder generically leads 
to magnetic frustration, because charge stripes act as 



antiphase domain walls for the AF order. Our results 
identify the hourglass as the excitation spectrum of the 
low-temperature cluster-spin-glass state which emerges 
due to the disorder-induced frustration. We discuss the 
ingredients required for this spin-glass physics and the 
implications for different material families. 

La^l^Srii^CoO/^. For definiteness, we choose to model 
the magnetism in La5/3Sri/3Co04, for which comprehen- 
sive INS data is available [9]. La2-a;Sra;Co04 shows sig- 
natures of stripe order for 0.3 < x < 0.6 [16, 17]. For 
X = 1/3 the Co^+ and Co^+ ions form a robust diagonal 
period-3 stripe pattern, with charge order setting in at 
Tco far above room temperature [17, 18]. Co^+ is in a 

= 3/2 high-spin state, while Co^+ is in a 5* = low- 
spin state [9, 19], and from the properties of the x = 
[20] and x = 0.5 [21] compounds it is known that the 
magnetic couplings - J between neighboring Co^+ and 
J' across a non-magnetic Co^+ - are both AF. For perfect 
charge order, this results in a stripe pattern as depicted in 
Fig. 1(a), exhibiting IC magnetism with magnetic Bragg 
peaks at Qs = (ttzLtt/S, 7r=F7r/3). [Stripes running along 
the opposite diagonal have Qs = (tt =b 7r/3, tt ± 7r/3).] 

Disordered charge stripes. The experiments in Ref. 9 
have proven substantial disorder effects to be present in 
La5/3Sri/3Co04. It is plausible to assume that those 
originate from static imperfections in the charge order. 
In contrast, dynamic fluctuations of the charge order are 
likely negligible for T <C Tco because of the robust in- 
sulating nature of the material which is re-enforced by 
spin-blockade eflPects [16, 17, 22, 23]. 

We therefore start by modelling stripes of Co^+ and 
Co^+ ions. We construct an Ising model for variables 
rii on the sites i of a square lattice, where = refers 
to a Co^+ and = 1 to a Co^+ ion. The model is 
chosen such that it has perfect stripe configurations with 
ordering wavevector Qc = (tt =b 7r/3, tt =F ^/3) [as in Fig. 
1(a)] and Qc = (tt ± 7r/3,7r ± 7r/3) as ground states at 
fixed filling (n) = 1/3, for details see Ref. 24. 

Monte-Carlo (MC) simulations of this model on L x L 
lattices (L < 48) are used to generate charge configura- 
tions away from this ordered state, i.e., with well-defined 
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Figure 1: (a) Perfect stripe order in the C0O2 planes of 
La5/3Sri/3Co04. Circles (arrows) represent non-magnetic 
Co^^ (magnetic Co^^) ions. Both magnetic couplings J 
and J' are AF. (b) Frustrated bonds caused by disorder in 
the charge sector, (c) Disordered charge configuration with 
5 and 2, obtained from the Ising model [24] for 

L = 48. Black (gray) squares correspond to Co^^ (Co^^) ions, 
(d) Charge structure factor Sc (q) for the same parameters as 
in (c) averaged over 10^ charge configurations. 

short-range order [24]. Representative results are in Fig. 
1(c) and Ref. 24, showing multiple stripe domains and 
various types of defects. We characterize these configu- 
rations by their correlation lengths (^^ ) parallel (per- 
pendicular) to the stripe direction, we obtain the from 
the charge structure factor Sc{q), Fig. 1(d), by a fit to 
a Lorentzian lineshape near Qc- We note that the pro- 
cedure of generating disordered stripes is not unique, as 
- for fixed ~ the types and distributions of defects 
will depend on details of both the Ising model and the 
MC protocol [24]. However, such details have remarkably 
little influence on the behavior of the spin sector. 

Spin model To model the magnetism of 
La5/3Sri/3Co04, we place localized 5* = 1/2 spins 
on the Co^+ sites of a 2D disordered stripe configuration 
as in Fig. l(b,c). The spins are assumed to interact via 
Heisenberg exchanges as in Fig. 1(a): 

ij a 

The first sum runs over the lattice sites with = ac- 
cording to a given charge configuration {n^}, and a = 
x^y^z. Guided by Refs. 9, 20, 21 we assume that the 
nearest-neighbor coupling J (which sets our energy scale) 
has the same form as in the undoped parent compound: 
= J (1 + e), jy = J, and = J (1 - d). The param- 
eters 6 and e control the spin anisotropy with 6 = 0.28 
and e = 0.013 [20, 21, 24]. Because the coupling across a 
Co^+ ion, J', is considerably smaller than J, we assume 



it to be isotropic with typically J' = 0.05 J. Interlayer 
couplings are neglected. 

Spin-glass ground state and excitations. For imperfect 
charge order, Hsp describes frustrated magnetism: As 
illustrated in Fig. 1(b), stripe defects generically lead to 
frustrated couplings between neighboring AF domains. 
This combination of disorder and frustration is expected 
to lead to spin-glass behavior, more precisely cluster spin- 
glass behavior, because the entities subject to frustrated 
couplings are not single spins, but domains of varying 
size. This spin glass displays short-range magnetic order 
(inherited from the perfect stripe state), with magnetic 
correlation lengths ^g- which we determine from the 
spin structure factor Ssiq)- In general, we find that 
follows [24] , but for a strongly disordered charge sector 
with < 4, Ss{q) is no longer peaked at the IC Qs, but 
instead displays a broad maximum at (tt, tt). This implies 
that the real-space period-3 structure is lost, and a spin 
glass with very short-ranged AF correlations emerges. 

To access the spin dynamics, we first determine locally 
stable classical states of Hsp for a given charge configura- 
tion {ui}: For classical spins 5*^, we perform MC simula- 
tions of Hsp utilizing single-site updates, combining the 
heat-bath and microcanonical (or over-relaxation) meth- 
ods. To efficiently sample all spin configurations, we also 
employ the parallel-tempering algorithm [25]. This pro- 
cedure is followed by a quench down to T = using 
a greedy conjugate-gradient algorithm. Our simulations 
confirm spin-glass behavior, as detailed in Ref. 24. 

We then calculate the excitation spectrum of Hsp us- 
ing linear spin- wave theory on finite lattices: Deviations 
from a classical state are represented by non-interacting 
bosons [24] modelling Gaussian magnetic fluctuations. 
The dynamic susceptibility is given by: 

X"(q,^)= [^^|(0|5"(<7)|^)|'<5(c.-c..)]^„, (2) 

a u 

where S'^{q) is the Fourier-transformed spin operator, |0) 
the magnon vacuum, and a single-magnon state with 
energy ujiy. We typically evaluate the average [• • • over 
80 spin states obtained from 40 charge configurations. 
Formally, such a calculation yields the zero-temperature 
excitation spectrum of a system with spatially disor- 
dered static stripes, but the results continue to apply 
at finite temperatures T and even for slowly fluctuating 
stripes, provided that both temperature and fluctuation 
frequency are small compared to the energies u of interest 

Charge disorder and hourglass spectrum. Fig. 2 shows 
our central result, namely the evolution of x"{^i^) with 
varying degree of stripe disorder for weak inter-stripe 
coupling (i.e. large J /J')- The perfect-stripe case 
Fig. 2(a) has sharp magnon peaks; at low energies gapped 
magnon branches (recall the magnetic anisotropy) dis- 
perse both inward and outward from Qs, and a sharply 
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Figure 2: Magnetic excitations, x"{^^^)^ of disordered stripes. From (a) to (e) the disorder in the charge sector increases, as 
characterized by the correlation lengths ^q, ^q. Each panel is divided at cj = cjs, with the lower [upper] portion showing data 
along q = {q^q) [{q^n)] in order to represent the high-intensity features. Parameters are J' = 0.05 J, L = 36, and (a) Perfect 
stripes, cjs = 0.88J5'; (b) 19, - 8, a;s = 0.82J6'; (c) ^ - 10, ^ 4, cjg = 0.70J5'; (d) ^ - 5, - 2, c^s = 0.50J5'; 

(e) 3, ^ 1, = 0.45J6' [24]. The S peaks in Eq. (2) were replaced by Lorentzians with width F = O.IJS. 



peaked large intensity maximum occurs at (7r,7r) corre- 
sponding to a saddle point of the spin-wave dispersion. 
In contrast, as disorder increases, Figs. 2(b,c), the ex- 
citations are broadened both in momentum and energy. 
The anisotropy gap is smeared, and the intensity at low 
energies is enhanced as compared to the clean case: Dis- 
order causes a strong redistribution of spectral weight 
from high to low energies [24]. Even more strikingly, for 
the strongly disordered case in Figs. 2(d,e), the outward 
dispersing branches at low energy have disappeared com- 
pletely. Instead, the weight is concentrated near (tt, tt) 
over a broad energy range, leading to an apparent "verti- 
cal" dispersion. Finally, in Fig. 2(e) the charge disorder 
is so strong that the IC signal at small u is hardly visible; 
here, Ss {q} shows a broad maximum at (7r,7r) [24]. 

A remarkable trend is that disorder causes - apart from 
broadening - a strong downward renormalization of the 
saddle-point energy: For otherwise identical parameters, 
ujs in Fig. 2(e) is half of that in the perfect-stripe case. 

Influence of model parameters. In the spin sector, the 
only parameters are the exchange constants J and J' 
and their anisotropics, whose values we fixed based on 
experimental data on La5/3Sri/3Co04 [9] and its parent 
compound [20]. Some uncertainty concerns the strength 
of the inter-stripe coupling, i.e., the ratio J/J^ In fact, 
this ratio varies considerably between different families 
of striped oxides, with values of 2 reported for nicke- 
lates (which do not show an hourglass dispersion) [26, 27] 
and 7-20 for cobaltates [9]. Spin- wave calculations in 
the disorder-free case [9, 26, 28] show that large J/ J' is 
crucial to the observation of an hourglass-like spectrum, 
whereas for small J/J^ the saddle point moves to the 
top of the spin-wave dispersion. We have performed cal- 
culations for disordered stripes with varying J/ J' which 
confirmed this analysis, for details see Ref. 24. 

Turning to the charge sector, our modelling here re- 
lies on a set of parameters and assumptions [24], with 
rather little experimental guidance available [18]. We 
have tested different schemes to introduce stripe disorder, 
in particular resulting in a varying degree of anisotropy 



^c/^C' We have found that this has rather little influ- 
ence on the magnetic behavior, such that our conclusions 
are independent of details of the stripe disorder [24] , and 
the only important parameter is ^c- 

Application to La2-xSrxCo04. A comparison of our 
results with the INS data of Ref. 9 indicates that charge 
disorder as in Fig. 2(d) is realized in La5/3Sri/3Co04. 
For these parameters, we show constant-energy cuts of 
x"{^i^) ill Fig. 3. We notice that the low-energy inelas- 
tic response {uo = 0.12JS) displays an anisotropic broad- 
ening - similar to the broadening of the Bragg signal. 
Fig. S5(a) [24], where ^/^^ ^ 2.3 - a behavior which 
is generic for our model. Significant low-energy intensity 
occurs also at (7r,7r), a feature that, in our calculations, 
is generically linked to an hourglass spectrum. 

For La5/3Sri/3Co04, we deduce JS ~ 25meV, for 
which we show an explicit comparison with the experi- 
mental data in Fig. 4. We find excellent agreement, apart 
from an intensity mismatch at 3 meV which is likely due 
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Figure 3: Constant-energy cuts of x"{^^^) with model pa- 
rameters as in Fig. 2(d): J' = 0.05 J, ^ - 5, - 2, L = 36. 
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Figure 4: Comparison of x\q^oj) for parameters as in 
Fig. 2(d) and JS = 25meV (solid) to the La5/3Sri/3Co04 
data (dots) taken from Fig. 3 of Ref. 9. Panels (a-c): 
q — (tt + ^,7r — ^); (d-f): q — (7r,7r + <J). Also shown are 
ad-hoc broadened perfect-stripe spectra (dashed) following 
Refs. 9, 24, with JS =18 meV adjusted to match the energy 
cjs of the saddle point. For both theories, the energy broad- 
ening was r = 1 meV, and an overall intensity factor and a 
constant background were adjusted by fitting the 14 meV data 
and then were used for all cj [29]. 



to the simplified treatment of spin anisotropies in our 
spin-1/2 model [24]. Fig. 4 also shows the result of ad- 
hoc momentum-space broadening of perfect-stripe spec- 
tra [24], as used in previous modelling of the data [9]. 
The most obvious failure of this approach concerns the 
missing intensity at energies below and slightly above 
the saddle point; it also overestimates the intensity of 
the upper spin-wave band which becomes visible near 
(tt, tt) at 35 meV. In addition, the mismatch of peak po- 
sitions at elevated energies. Fig. 4(d,e), illustrates that 
the real hourglass has approximately linearly dispersing 
branches above the saddle point, while the perfect-stripe 
case yields a parabolic dispersion - a difference repro- 
duced by our theory which becomes clear upon compar- 
ing panels (a) and (d) of Fig. 2. 

We conclude that the charge stripes in La5/3Sri/3Co04 
are rather disordered, with an estimated correlation 
length oi ~ 5; this explains why charge order has not 
been observed in scattering experiments. Our analysis 
predicts the magnetism in La5/3Sri/3Co04 to be that of 
a cluster spin glass. Indeed, initial signatures for glassy 
behavior have been detected [30], and a more thorough 
characterization is called for. 

Finally, we propose that disordered stripes are also 
present in La2-ccSra^Co04 at x < 0.3: Smaller x im- 
plies larger stripe spacing, with increased sensitivity to 
quenched disorder. Indeed, both a broad Bragg signal 
[17] and broad low-energy excitations [31] have been de- 
tected in Lai.7Sro.3Co04 around (7r,7r), consistent with 
the magnetism of strongly disordered stripes. 

Other transition- metal oxides. For overdoped single- 
layer manganites with ordering wavevector Qs = 



(27r ± ^, TT ± ^) with S = O.IGGtt, INS measurements have 
recently detected hourglass spectra, in particular at ele- 
vated temperatures (while low-T spectra displayed out- 
ward dispersing branches) [10]. The hourglass has been 
linked to a short magnetic correlation length in this 
regime. Our model calculations qualitatively confirm this 
interpretation: For intermediate J/J^ the coupling be- 
tween AF domains is robust at low T, but weakened at 
elevated T, such that the advertised hourglass excitation 
spectrum of weakly correlated AF clusters is observed. 

Finally, it is interesting to link our study to cuprate su- 
perconductors. In fact, the physics of disordered stripes, 
either slowly fluctuating or disorder-pinned, has been in- 
voked to explain a variety of experimental observations 
including INS data [14, 15]. The hourglass spectra ob- 
served in cuprates show a clear absence of outward dis- 
persing low-energy branches, and moreover display a pro- 
nounced piece of vertical dispersion in the so-called pseu- 
dogap regime [2, 3, 6, 15]. Assuming that our modelling 
qualitatively applies to cuprates as well, we can link these 
experimental results to strongly disordered stripes. We 
also note that very underdoped (but still superconduct- 
ing) cuprates have been found to display a quasieleas- 
tic "central mode" at (tt, tt) which has been interpreted 
in terms of glassy short-range order [32]. Considering 
that decreasing the doping level in cuprates below 1/8 
weakens stripe order (such that it is more susceptible to 
quenched disorder), this nicely matches our results: For 
strong stripe disorder, Ss{q) looses its IC structure and 
instead displays a broad peak at (7r,7r) [24]. 

Summary. Our detailed modelling of disordered stripes 
in La5/3Sri/3Co04 reveals a link between the celebrated 
hourglass excitation spectrum and a cluster-spin-glass 
magnetic ground state, caused by frustrated couplings 
between locally ordered AF domains. These phenomena 
appear common to a broad class of doped oxides. 
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I. GENERATING CHARGE CONFIGURATIONS 

The charge order of Co^+ and Co^+ ions governs the 
magnetism of La2-a;Sra;Co04. As this charge order sets 
in above room temperature, we assume the charge order 
to be static on the time scale of magnetic fluctuations, 
such that the spin sector can be analyzed - in the spirit of 
a Born-Oppenheimer approximation - for a set of frozen 
charge conflgurations. Those conflgurations correspond 
to imperfect stripes, with imperfections likely caused by 
crystalline defects. 

We obtain such charge conflgurations from numerical 
simulations of an appropriate Ising model. Consider vari- 
ables Ui on the sites z of a square lattice, where = 
refers to a Co^+ ion, while = refers to a Co^+ ion. 
The rii are subject to a Hamiltonian of the form 



He = Y.^^Y. 



(SI) 



with a particle density of = adequate to de- 

scribe the compound La5/3Sri/3Co04, where N = L x L 
is the number of sites, and periodic boundary conditions 
are applied. The Vm represent density-density interac- 
tions as depicted in Fig. SI, which are chosen such that 
the ground state of He is given by the diagonal stripe 
state established for La5/3Sri/3Co04, Fig. SI. Specifl- 
cally, Vi is a nearest-neighbor repulsion, V2 < a second- 
neighbor attraction, and V3 > a third-neighbor repul- 
sion, while V4^5 are longer-range interaction used to mod- 
ify the defect structure (see below, for most of the calcu- 
lations V4,5 = 0). 

We generate short-range-ordered charge conflgurations 
using Monte-Carlo (MC) simulations of the Ising model 
in Eq. (SI). As noted above, the relevant charge disorder 
is not generated by thermal fluctuations at temperatures 
near or above a putative charge-ordering transition (be- 
cause this temperature is far above room temperature in 
La2-a:Sra;Co04), but instead arises from structural de- 
fects (e.g. the dopant distribution). As details of this 
quenched disorder are not known, and in order to keep 
the set of input parameters minimal, we choose to em- 
ploy the model Eq. (SI) without quenched disorder, but 
instead extract conflgurations away from thermal equi- 
librium. 

To this end, the variables {n^} are initialized by ran- 
domly distributing N/3 charges on N sites. One MC step 
consists of randomly selecting two sites and exchanging 
their charges, accepting or rejecting this move according 



(a) (b) 

oo«oo« o#oo#o 
00 •oo#oo 

)0#00# o«oo«o 

^ foo^ofo •oo#oo 

)0#00i 0«00«0 

oojf 00 •oo^oo 

V 



Figure SI: Perfect diagonal period-3 stripe patterns in the 
C0O2 planes of La5/3Sri/3Co04, where black and gray circles 
correspond to n^ = 1 (Co^^) and n^ = (Co^^), respectively. 
These configurations are ground states of the Ising model in 
Eq. (SI), with density-density interactions Vi > 0, 1^2 < and 
V3 > as shown in panel a). The longer-range interactions 
¥4,5 will be utilized to modify the structure of charge-order 
defects away from the ground state. 



to the rules of the standard Metropolis algorithm at a 
finite effective temperature T. We work at a low tem- 
perature corresponding to the charge-ordered phase and 
select configurations {n^} after a fixed number Nuc of 
MC steps. Thus, for large A^mc the system reaches a 
state with long-range order, while smaller A^mc result in 
short-range-ordered states. We assume that these states 
sensibly represent the situation in La2-a^Sra;Cu04, be- 
cause working at low T ensures to have well-ordered re- 
gions separated by localized defect structures, as can be 
expected for a material far below its charge-ordering tem- 
perature. Observables are finally obtained by averaging 
over a large number of charge configurations, all obtained 
from runs with different initial conditions - this mimics 
the average over a large disordered sample. 

To quantify the charge disorder, we calculate the 
charge correlation length ^c, which we obtain from the 
static charge structure factor 



= ^[J2{l-n,){l-nj)e- 



, (S2) 



where Ns = 2N/3 is the number of Co^+ ions (i.e. spins), 
and denotes the average over charge configura- 

tions. For short-range order, Sc {o} is typically peaked at 
four wavevectors Qc (see Fig. 1(d) of the main paper and 
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Figure S2: Evolution of the charge correlation length 
with A^MC, i.e., along the MC simulation of the Ising model 
Eq. (SI) for two different system sizes L. The parameters 
= 4, 1^2 = -2, 4, 14,5 = and T = 1 correspond 

to the charge-ordered phase, such that reaches the system 
size for large Nmc- The inset shows the evolution of the 
energy E = [Hc]av along the MC simulation, normalized to 
the energy of the perfectly ordered state Eq. 



Figure S3: Spin correlation lengths, and of the clas- 
sical ground state of Tisp, Eq. (S3), as function of the charge 
correlation length which characterized the stripe disorder. 
Clearly, ^| ~ while is smaller by a factor of 2-2.5. 



II. SPIN MODEL 



Figs. S4, S6): the peaks at Qc = (27r/3, 27r/3) and Qc = 
(47r/3, 47r/3) correspond to stripe domains like those 
in Fig. Sl(a), while peaks at Qc = (27r/3, 47r/3) and 
Qc = (47r/3, 27r/3) belong to stripes running along the 
opposite diagonal, Fig. Sl(b). The width of the peaks 
is anisotropic, thus we distinguish correlation lengths 
and parallel and perpendicular to the stripe direction, 
respectively. The evolution of with the number of 
steps, A^MC, within our non-equilibrium MC procedure 
is shown in Fig. S2. As anticipated, the ordered state 
(with ^ L) is reached for large TVmc, while smaller 
values, A^MC = lO'^-lO^, are suitable for obtaining charge 
configurations with ^ 2-10 lattice spacings. 

Individual charge configuration from our MC simula- 
tions are shown in the top rows of Figs. S4(b)-(e) and 
S6. We see that the protocol generates finite-size stripe 
domains of both orientations, as intended. For a large 
range of parameters with V4,5 = as in Fig. S4 we gener- 

ically find > ^q, i.e., the diagonal arrangement of 
neighboring charges tends to be more stable than the 
spacing between the diagonal charge lines. Sq {q) also 
reveals that for very large disorder. Fig. S4(e), signifi- 
cant weight appears away from Qc- This effect is even 
more pronounced in the spin sector, as will be discussed 
below. 

Charge configurations with < can be generated 
on purpose, e.g., by choosing the longer-range couplings 
^4,5 / 0, which stabilize the stripe spacing and cut charge 
lines, respectively. Results are shown in Fig. S6, which 
now show a visibly different defect structure in the charge 
sector. Remarkably, the behavior on the spin sector is 
rather insensitive to this change (see below). 



As explained in the main text, the magnetism of 
La5/3Sri/3Co04 is assumed to be well described by im- 
mobile magnetic moments that are located on the Co^+ 
sites (i.e. those with = 0) of a disordered 2D stripe 
configuration. 

The Co^+ ions are in a 3d^ high-spin configuration with 
5* = 3/2. A full modelling of these high-spin ions and 
their interactions has been put forward in Refs. 1,2 - this 
involves spin and orbital angular momentum, crystal- 
field and spin-orbit effects, as well as (nearly isotropic) 
exchange interactions. As the local Hilbert space in 
this model contains 28 states, it is more efficient to em- 
ploy a reduced effective model instead: As shown in 
Refs. 1,2, the low-energy part of the excitation spectrum 
of La2-a^Sra^Cu04 both at x = and x = 0.5 can be well 
captured by a 2D spin-1/2 spin-only model. 

Consequently, we adopt such a model for our purpose, 

i,j a 

and choose the nearest-neighbor J as in Refs. 1,2 to 
be = J(l + e), jy = J, and = J{l-S), with 
S = 0.28 and e = 0.013. This exchange anisotropy mim- 
ics the effect of the original single-ion anisotropy, result- 
ing from the combination of crystal-field and spin-orbit 
interactions. It is chosen such that the ordered AF state 
has spins in the x direction, corresponding to the crystal- 
lographic a axis, as in experiment. Following Ref. 3 we 
neglect the anisotropy of the smaller coupling J' across 
a non-magnetic Co^+ ion. 

As discussed in the main text, for imperfect charge or- 
der, the classical ground state of Hsp is a cluster spin 
glass with well-defined short-range magnetic order, see 
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also Sec. VI below. In Fig. S3, we demonstrate the cor- 
relation between and ^c- In particular, we obtain 
^ system sizes considered. 

III. SPIN- WAVE CALCULATION 

We calculate the excitation spectrum of Hsp (which in- 
volves spatial disorder) using real-space linear spin- wave 
theory^^ on finite L x L lattices with periodic boundary 
conditions. This is done by expanding around a locally 
stable classical state (or saddle point) of Hsp as deter- 
mined by the spin MC procedure described in the main 
text. The relevant classical states are in general non- 
collinear, with the spin directions of the Si characterized 
by angles Oi and We start by introducing a new ref- 
erence frame in which points locally in the direction 
of the magnetization 




cos^^cos^^ 
cos^isin^^ 
—s'mOi 



—sinq 

cos^. 




sin^^cosi 
sin^^sin- 
cos^- 




where uji, is the eigenfrequency associated to the uth 
mode, and 
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(S8) 



with = (ai,a2,--- ,aArs), = (&i,^2,-" ^^a^s)- U 
and V are A^s x matrices describing the transforma- 
tion to this diagonal basis. ^'^ 

Finally, we evaluate the rotationally averaged dynamic 
susceptibility, X^\q^oj)^ in the standard one-magnon ap- 
proximation, see Eq. (2) of the main paper. For a single 
realization of disorder and using the spin-wave transfor- 
mation as described above we have 

X"{q,uj) = ^^|M,"(5-)|2<5(^-^,). (S9) 

V a 

Here a = z and (q) represents the matrix ele- 
ment (0|5'" (q) \u) of the Fourier-transformed spin oper- 
ator with the eigenmode \iy): 



L_^e-«"'-«M«(r-;), (SIO) 



In this rotated frame S\ we apply the linearized Holstein- 
Primakoff transformation: S^^ = \/2Siai, S^' = 
\/2Sia\^ and = Si — ajai, where 5^ = 1/2 is the spin 
size at site i. The bilinear piece 7^2 of the Hamiltonian 
then reads 



where 



with the coefficients 



Myifi) = W y [t/i^ (cOSl9iCOS^ 

+ V^^ {cos6iCos(j)i — isin^ 



- isin^^) + 



^2 = E ^^Maj + \ E [BiAa] + h.c) , (S4) (n) = ^ ^ [C/,. (cos^.sin, 



icosd. 



i{F^y-F^y)], 



Bi 



(S5) 

{F^y + F^y)]. (S6) 



+ V*, (cos^isin(/>j + icos(j)i)\ 
2 



(Sll) 

(S12) 
(S13) 



As momentum is not a good quantum number in a disor- 
dered system, the matrix element (g), for any given 
g, will in general be non-zero for all eigenmodes v. 



Here, the functions F^-'^ read 



IV. CHARGE DISORDER AND HOURGLASS 
SPECTRUM: ADDITIONAL RESULTS 



rnxx 



= JfjCos6iCos6jCos(t)iCos(j)j + 

+ jyjCOs6iCOs6jS'm(t)iSm(t)j + J^^sin^^sin^j, 

= J^^sin^isin^j + JfjCOS(j)iCos(j)j, 



Jfj ^hiOi ^hiOj cos(/)^cos(/)j 

jy 



+ J|:sin6>^sin6>jsin(/)isin(/)j + J?cos^^cos^j, 



JfjCosOiCos(j)isin' 



Jfj cosOi s'm(j)i cos(j)j . 



The bilinear Hamiltonian Eq. (S4) can be nu- 
merically diagonalized using a bosonic Bogoliubov 
transformation,^ and the resulting Hamiltonian reads 
(up to a constant) 



no 



(S7) 



The evolution of the spin excitations with increasing 
charge disorder is shown in Fig. 2 of the main paper. 
Here we complement this by Fig. S4 which displays snap- 
shots of the charge configuration as well as Sc{q) and 
5's(^ for the same parameters as in Fig. 2. The figure 
reveals a very similar progressive and anisotropic broad- 
ening in Sc{q) and Ss{q) for weak and moderate disorder. 
This is different for strong disorder. Fig. S4(e), where 
the IC signal in Ss {q) has almost disappeared in favor of 
a broad peak at the AF wavevector (7r,7r), while Sc{q) 
still displays a clear IC structure. This corresponds to 
a state where disorder-induced frustration has destroyed 
the period-3 spin correlations, leaving behind weakly cor- 
related local AF clusters - this is further illustrated in 
Fig. SIO below. 
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Figure S4: Numerical results for individual charge configurations rii (top), the charge structure factor Sc{q) (middle), and 
the spin structure factor 5's(^ (bottom) for disordered stripes with L = 36. From (a) to (e), the degree of charge disorder 
increases as in Fig. 2 of the main paper. For the charge sector, the parameters are Vi = 4, V2 = —2, V3 = 2, V4,5 = 0, T = 1, 
and (a) perfect stripes, (b) A^mc = 2 x 10^ 19, - 8, (c) A^mc = 0.5 x 10^ 10, - 4, (d) A^mc = 0.1 x 10^ 

~ 5, — 2, (e) A^MC = 0.025 x 10^, ~ 3, ^ 1- The spin excitation spectra corresponding to the situations in panels 
(a) to (e) are displayed in Fig. 2 of the main paper, which also specifies the parameters of the spin model 1-Lsp- 



A. Spin-sector parameters 
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We have performed calculations for disordered stripes 
with varying J/J^ with sample results in Fig. S5. These 
confirm that J/J^ decisively influences the form of the 
magnetic spectrum: While an hourglass emerges for 
J/ J' = 20, for there is neither a large-intensity 

feature at (tt, tt) nor a "vertical" dispersion. Instead, large 
weight is concentrated at low energies close to Qs, even 
though the outward dispersing branches are weakened 
by disorder. Moreover, the structure factor Ss{q) reveals 
that the spin sector is less sensitive to charge disorder for 
smaller J/J^: stronger inter-stripe coupling stabilizes the 
IC state against glassy behavior. 

We note that J /J' is often deduced via the experi- 
mentally measured location of the saddle point, uJs/{JS). 
Our results in Fig. 2 of the main paper show that care is 
required, because uJs/{JS) also depends on the degree of 
stripe disorder. 
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Figure S5: Influence of the coupling ratio J/ J' on the mag- 
netism, for charge-stripe parameters as in Fig. 2(d) of the 
main paper. (a,b) Ss{q) and (c,d) low-energy part of x"{q^^) 
along q= {q,q) for (a,c) f = 0.05J and (b,d) J' = 0.15J. 



B. Charge-sector parameters 

As mentioned in Sec. I, the disordered stripes obtained 
from our MC procedure quite generically display > 
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Figure S6: Numerical results for disordered stripes as in 
Fig. S4, but here for stripe parameters with ^f. < «fc • From 
top to bottom: charge configuration n^, charge structure fac- 
tor Sc{q), spin structure factor Ss{q), and the low-energy part 
of the magnetic excitations, x'{Qj^): along q = {q,q)- The 
parameters for the charge sector are L = 36, Vi = 2, V2 = 0, 
V3 = 2, ^4 = -2,1/5 = 0.5, T = 1, and (a) A^mc = 4 x 10^, 
^ 5, - 10, (b) iVMC = 0.5 X 10^ ^ 2, - S.^For 
the spin sector we have J/ J' = 20 and the S peaks in x'\q^ ^) 
were replaced by Lorentzians with width F = O.IJS*. Remark- 
ably, Ss{q) displays anisotropic broadening with ^| > ^s"? Re- 
spite < <Jc in the charge sector. Consequently, the spin 
excitations here are qualitatively similar to that in Fig. 2 of 
the main paper (for a comparable degree of charge disorder). 



- this applies to all cases in Fig. S4. Results for the op- 
posite situation, < , are shown in Fig. S6 - those 
were obtained using V^^^ 7^ in He, Eq. (SI). Remark- 
ably, even for this kind of disorder configuration we find 
^1 is which thus appears universal. This can be ra- 
tionalized considering that we work in the limit J ^ J': 
Here, intact AF domains will show nearly perfect order. 



while correlations between AF domains, mediated by the 
small J', can easily be destroyed by spatial disorder. As 
the geometry implies that J stabilizes correlations along 
the stripe direction, this explains ^| > ^g". It is then 
not surprising that the evolution of the spin excitation 
spectrum in the case < ' shown in the bottom row 
of Fig. S6, in qualitatively similar to that in the case 
^C' Eig. 2 of the main paper: While both high- 
intensity features at Qs and outward dispersing branches 
at low energy are still visible for weak disorder, those 
tend to be suppressed for strong disorder, where instead 
an apparent "vertical" dispersion appears at intermediate 
energies. 

An interesting feature in Fig. S6(b) is a considerable 
enhancement of the intensity at (tt, tt) even for ener- 
gies well below the saddle point, coexisting with low- 
energy weight near Qs- This can be traced back to the 
momentum-space structure of 5's (0, which shows peaks 
both at Qs and (7r,7r), arising from the coexistence of 
relatively well-defined period-3 structures and small AF 
clusters. (We have not observed this pronounced coexis- 
tence for disordered stripes with > ' implying that 
microscopic details of the charge sector play some role 
here.) A further increase of the charge disorder tends 
to transfer spectral weight from Qs to (7r,7r), similar to 
Fig. 2(e) of the main paper. 



V. COMPARISON TO BROADENED 
PERFECT-STRIPE SPECTRA 

In order to further illustrate the importance of a proper 
treatment of stripe disorder, we provide a comparison of 
our results to those from spin-wave calculations for per- 
fectly ordered stripes, where disorder is simply accounted 
for by an ad-hoc broadening in momentum space. This 
procedure has been employed previously in the literature: 
For instance, in Ref. 3 reasonable agreement with the ex- 
perimental data concerning the intensity distribution in 
momentum space was achieved. 

We have implemented the ad-hoc broadening for the 
spin-1/2 model in Eq. (S3) following Ref. 3, by convolv- 
ing x"{Q->^) at a fixed energy uj with a two-dimensional 
anisotropic Gaussian. This introduces two fit parameters 
- the momentum-space widths 7|| and 7^ - which mim- 
ick stripe disorder and thus replace the physical parame- 
ter of our calculation. Results for broadened perfect- 
stripe spectra with J / J' = 20 are shown in Fig. S7, where 
7|| and 7^ were chosen as to match the momentum width 
of the experimentally observed excitations at low energy 
(3meV in Fig. 4 of Ref. 3). (The same parameters are 
used for the dashed curves in Fig. 4 of the main text.) 
As a result, the momentum-space intensity distributions 
are roughly similar to that obtained from the full calcu- 
lation displayed in Fig. 3 of the main paper (and also 
roughly similar to that of the experimental data in Fig. 4 
of Ref. 3). 
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Figure S7: Constant-energy cuts of x"{Qi^) ^ls a function 
of momentum q as in Fig. 3 of the main paper, but now for 
perfect stripes with J / J' — 20 and an anisotropic ad-hoc 
broadening in momentum space of 7^ = O.IStt and 7|| = 
O.OOtt. These results are broadly similar to that in Fig. 4e-h 
of Ref. 3, where, however, the full crystal- field model for the 
Co^^ ions (instead of an effective spin- 1/2 model) was used. 




CO/JS 

Figure S8: Momentum-integrated spectral weight as function 
of energy cj, comparing the result from a full simulation of 
disordered stripes (with 5, — 2 as in Fig. 2d of the 

main paper) to that of perfect stripes (additional momentum- 
space broadening does not influence this curve) . In both cases 
J I J' = 20, and a Lorentzian energy broadening of F = O.IJS* 
has been employed. The arrows denote the energy position 
cjs of the (7r,7r) saddle point, and the intensities have been 
rescaled to match at cjs- The intensity peaks in the range 
1.7 — 2. US originate primarily from the upper of the two 
spin- wave branches.^ 



However, a closer inspection at the energy depen- 
dence of the intensities reveals strong differences. This 
is already seen in the comparison at selected energies 
as shown in Fig. 4 of the main paper, and becomes 
clear when considering the momentum-integrated inten- 
sity x"{'^)i Fig. S8. While the perfect-stripe-based re- 
sult has rather little intensity at low energies and a 
sharp saddle-point feature at uos^ the result for disor- 
dered stripes has large low-energy weight (note that even 
the anisotropy gap is almost completely smeared) and a 
broad saddle-point feature. As a result, the latter calcu- 
lations yields only moderate intensity variations between 
0.1 75 and 2.5 JS^ in accord with experiment (compare 
Fig. 4 of the main paper). In contrast, much larger inten- 
sity variations with energy arise from perfect stripes. In 
addition. Fig. S8 underlines the downward renormaliza- 
tion of ujs and also the overall increase of the bandwidth 
caused by disorder. Therefore, extracting JS (or J/ J') 
by matching theoretical and experimental features of the 
spectra (like the saddle-point location us) requires care. 



VI. SPIN-GLASS SIGNATURES 

We can characterize spin-glass behavior in our MC sim- 
ulations for the spin sector by monitoring the specific 
heat C and the magnetic correlation length as func- 
tion of temperature T, both shown in Fig. S9. Panel (a) 
shows that the specific-heat peak, which occurs near the 
mean- field Neel temperature in the perfect-stripe case, 
is strongly broadened and shifted to a lower Tf with in- 
creasing stripe disorder. This effect is more pronounced 



for weakly coupled stripes, i.e., large J/J^ As can be 
seen in panel (b), the correlation length grows upon 
cooling, in particular in the temperature regime near Tf , 
but remains much smaller than the system size as T ^ 0. 
Again, this is more pronounced for large J/J^ which 
also causes a larger anisotropy as measured by the ra- 

tio 4/^^. 

In addition, we have also calculated the spin-glass cor- 
relation length ^sG from the correlations of the Edwards- 
Anderson order parameter.^' This correlation length 
strongly increases below Tf as well and reaches values of 
the order of the system size as T ^ (not shown). Im- 
portantly, ^sG is always larger than (except for very 
large T). 

All these results are consistent with the behavior of a 
spin glass with short-range magnetic order and a mean- 
field freezing temperature Tf. Clearly, glassiness in- 
creases with increasing charge disorder and with increas- 
ing J/J^ as indicated by a broader freezing regime and 
smaller ^s- For large J/J^ the glassy state consists of 
ordered AF clusters (formed by intact stripe segments of 
Co^-h ions) which are weakly correlated among each other 
(see also Fig. SIO) - this justifies to call this a cluster spin 
glass. Note that the ground state is expected to remain 
glassy even in the limit of J ~ J', where the strongly 
disordered system corresponds to a vacancy-doped frus- 
trated spiral magnet which displays more conventional 
spin-glass (instead of cluster-spin-glass) behavior. 

We note that, in our two-dimensional spin model Hgp, 
a true spin-glass transition does not exist. Therefore, 
we refrain from a detailed finite-size scaling analysis of 
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the thermodynamic and magnetic properties'^' of Hsp - 
such an analysis would also be complicated by the fact 
that we have no straightforward way to introduce pre- 
cisely the same amount of microscopic disorder at dif- 
ferent L due to the use of a non-equilibrium MC pro- 
tocol. We have, however, investigated the L depen- 
dence of our results for the spin structure factor Ss{q) 
at fixed (which is indeed a very good estimator of 
disorder). For Ss{q) we essentially find no dependence 
on the system size in the range L = 24-48. In particu- 
lar, Ss{Qs)/Ns ^ as l/7Vs ^ (see the inset of Fig. 
S9(b)) indicating once more static short-range order with 
a vanishing magnetic order parameter. Similarly, the sus- 
ceptibility X^\q,uj) displays no signs of size dependence 
(not shown). 



T/JS 



Figure S9: Characterization of glassy behavior in the spin 
sector, as obtained from the MC simulations of Tisp for 
L = 36. (a) Specific heat C as function of temperature, for 
different J /J' (see legend) and different amounts of charge 
disorder: circles: 5, <Jc ~ 2, squares: ^ 10, <fc ~ 5, 

solid line: perfect stripes, (b) Spin correlation lengths 
as function of temperature for different J/ J' and ^ 5, 
^ 2. Here, the results at T = have been obtained from 
locally stable classical states. Error bars are typically smaller 
than the symbol size. Note that, for higher T (where no data 
points are shown), Ss{q) does not show well-defined peaks 
at Qs but instead a broad maximum at (7r,7r), i.e., we ob- 
serve a temperature-dependent onset of incommensurability. 
Inset: T = spin structure factor Ss{q) at one of the ordering 
wavevectors Qs normalized by the number of spins A^s as a 
function of 1/Ns. For both curves we have J /J' = 20 and 
different amounts of charge disorder: circles: ^(|?^5, ^c~2, 
squares: ~ 10, ~ ^^ 



Finally, we illustrate the cluster-spin-glass state by 
showing specific saddle points of Hsp in Fig. SIO: these 
are portraits of the spin configurations of locally stable 
classical states. The parameters correspond to moder- 
ately and strongly disordered stripes as in Fig. 2(d,e) 
and Fig. S4(d,e). The configurations nicely show clus- 
ters with local AF order and essentially random relative 
orientations. In particular, the states are non-collinear, 
i.e., the small easy-axis anisotropy e is less effective in 
the presence of disorder. 
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Figure SIO: Real-space spin configurations in classical ground states of Tisp with L = 36, for the same parameters and charge 
configurations as in Fig. S4(d) (^^ ^ 5, — 2 (a)) and Fig. S4(e) (^^ — 3, < 1 (b)). The dots show the non-magnetic 
sites, while the arrows show the x and y components of the Si (the z components are tiny due to the strong anisotropy caused 
by S). 
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